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ABSTRACT 

We present the mid-infrared (MIR) spectra obtained with the Spitzer InfraRed Spectro- 
graph (IRS) for a sample of 52 sources, selected as infrared luminous, z > 1 candidates in 
the Extragalactic First Look Survey (XFLS). The sample selection criteria are fiA^m <^ 0.9mJy, 
uf u (2A^m)/iyf u (8^m) > 3.16 and vf l/ (2£fj,m)/vf„(0.7fj,m) > 10. Of the 52 spectra, 47 (90%) 
produced measurable redshifts based solely on the mid-IR spectral features, with the majority 
(35/47 = 74%) at 1.5 < z < 3.2. Keck spectroscopy of a sub-sample (17/47) agrees with the 
mid-IR redshift measurements. The observed spectra fall crudely into three categories — (1) 
33% (17/52) have strong PAH emission, and are probably powered by star formation with total 
IR luminosity roughly a factor of 5 higher than the local starburst ULIRGs. (2) 33% (17/52) 
have only deep silicate absorption at 9.8 fim, indicative of deeply embedded dusty systems. The 
energetic nature of the heating sources in these systems can not be determined by these data 
alone. (3) The remainder 34% are mid-IR continuum dominated systems with either weak PAH 
emission and/or silicate absorption. This third of the sample are probably AGNs. We derived 
monochromatic, rest-frame 5.8 /im, continuum luminosities (yL v ), ranging from 10 10 ' 3 — 10 12 ' 6 Lq. 
Our spectra have MIR slope a^^i^tmi <L 2.1, much redder than the median value of 1.3 for the 
optically selected, Palomar-Green (PG) quasars. From the silicate absorption feature, we esti- 
mate that roughly two-thirds of the sample have optical depth rg.8 Mm > 1. Their and 
.Lir suggest that our sample is among the most luminous and most dust enshrouded systems of 
its epoch. Our study has revealed a significant population of dust enshrouded galaxies at z ~ 2, 
whose enormous energy output, comparable to that of quasars, is generated by AGN as well as 
starburst. This IR luminous population has very little overlap with sub-mm and UV-selected 
populations. 
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1. Introduction 

Ultra-luminous Infrared Galaxies (ULIRGs, 
Lg-iooo/im > 1 12 L(T)), while relatively rare in the 
local universe ( Soifer et al. Il987t ISanders fc Mirabell 
1996), are far more common at high redshifts, 
as suggested by observations from ISO , (sub)- 
millimeter, and more recent l y Spitzer (|Elbaz et al 



put from ULIRGs makes them sig nificant contrib- 
utors to global luminosi ty density (jGuiderdoni et al 



19981 : 
found 
tions 



Dole et al. 20061 ). Detailed studies have 
evidence for the evolutionary connec- 
between ULIRGs and the formation of 



quasars and elliptical galaxies (iDasvra et al. 



2006tlVeilleux et al. l200a I Alexander et al.ll2005l ). 



19991: IChapman et al.ll2003t iLe Floc'h et alj |2005: 



Daddi et al.ll2005h . The tremendous energy out 



The unprecedent ed sensitivity of th e InfraRed 
Spectrograph (IRS: iHouck et ah 2004 ) on Spitzer 
has made it possible for the first time to observe 
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galaxies at z > 0.5— 3, a virtually unexplored ter- 
ritory for mid-IR spectroscopy. The rest-frame 
3 — 25 fim wavelength region is home to a rich 
suite of spectral features, essential for understand- 
ing the physical properties of ULIRGs. Low res- 
olution, mid-IR (3-40/xm) spectra of Luminous 
Infrared Galaxies (LIRGs) and ULIRGs consist of 
broad Aromatic emission (3.3, 6.2, 7.7, 8.6, 11.2, 
12.7 /im) and silicate absorption at 9.7 /im plus 
various strong ionic or molecular lines, such as 
[Aril], [Nell], [Ne lll], H 2 , CO etc., superposed on 
a red continuum drlouck et al.l[2004t Armus et al 



120041 : iBrandl et al.ll2004l : ISnoon et al.ll2004h . The 

Aromatic features are usually ascribed to Poly- 
cyclic Aromatic Hydrocarbons (PAH) , with < 200 
atoms and fluctuating effective te mperatures re- 
spon ding to stochastic heating (jPuget fc Leger 
1989). These PAH features are prominent in star- 
forming systems, but are reduced and modified in 
high-intensity starbursts, eventually disappearing 
in AGN systems. They provide both an indication 
of the energy source heating the dust and red- 
shift estimates for sources tha t are completely ob- 
scured at shorter wavelengths dGenzel fc Cesarskv 



2000; iDrainefeOoilLaurent et al.ll2000l: iTran et al 



20011: iRigopoulou et alj 11999 : IVoitl 119921 ). The 



mid-IR continuum is associated with very small 
dust grains, which are more resilient to photo- 
destruction than the PAHs (Draine & Li 2001). 
The mid to far-IR continuum ratio constrains the 
amounts of hot and cold dust. Finally, the depth 
of silicate absorption is a quantitative indicator of 
the mid-IR dust opacity along the line-of-sight. 

In this paper, we report the observational re- 
sults from the Spitzer GO-1 program (PID: 3748), 
which provides an initial characterization of mid- 
IR (7-38/xm), low resolution spectroscopic prop- 
erties of a sample of high-redshift, infrared lumi- 
nous galaxies selected in the Extragalactic First 
Look Survey (XFLS). This program aims to ad- 
dress the basic spectral properties of high-z, IR 
luminous galaxies. This paper is the first in a se- 
ries of papers from this survey. Its main scope is 
to present the target selection, observations (§2), 
spectral reduction method (§3), observed spectra 
and redshift distribution (§4). We discuss the im- 
plication of our findings in §5. The detailed anal- 
yses of the spectral properties — PAH strength, 
mid-IR opacities and continuum slopes — and the 
cosmological implication for the high redshift in- 



frared luminous populations are presented in two 
subsequent papers (Sajina et al. 2006a, PAPER 
II hereafter; Sajina et al. 2006b, in preparation). 

Throughout the paper, we adopt the standard 
O M = 0.27, O\ = 0.73 and H = 71 kms^ 1 Mpc" 1 
cosmology. 

2. Target Selection and Observations 
2.1. The Target Sample 

We selected a total of 52 targets in the Spitzer 
XFLS3 over an area of 3.7deg 2 , according to 
the following criteria: (1) > 0.9mJy; (2) 

#(24,8) = log l0 (uf v (24fj,m)/vf u (8fj,m) > 0.5; (3) 
#(24,0.7) = log 10 ( i // 1 ,(24Mm)/^/ l/ (0.7/im) > 1.0. 
We used the Spitzer and #-band catalogs pub- 
lished in Lacy et al. (2005), Fadda et al. (2006), 
and Fadda et al. (2004). The 24 ^m flux density of 
0.9mJy cut off was chosen so that moderate signal- 
to-noise ratio (S/N) spectra could be obtained in 
a reasonable amount of tel escope time. Using th e 
infrared spectral template ( Chary fc Elbaz 112001 ). 
the 0.9 mJy at 24 {im limit crudely translates to 
Lib. = ^l-iooopm ~ 10 12 L Q and 1O 13 L at z = 1 
and 2 respectively. As discussed in detail in §4.4, 
the high luminosity of our sources determines the 
spectral properties observed among our sample. 

The 24-to-8/im and 24-to-R color cuts are to 
select galaxies with strong PAH or steep mid- 
IR continuum at z > 1. This color selection 
technique has been discussed in detail by Yan 
et al.(2004), and more recently, by Brand et al. 
(2006). The i?(24,0.7) > 1 color cut corre- 
sponds to rest-frame z// l/ (8/im)/i// iy (2333A) > 10 
at z ~ 2, selecting infrared bright and optically 
faint sources. At f v (24(j,m) = lmJy, #(24,0.7) > 
1 implies # magnitude fainter than 22.6 magni- 
tude (Vega). The 24-to-8/xm color criterion picks 
out sources with strong PAH and/or steeply rising 
mid-IR continuum. #(24, 8) > 0.5 corresponds to 
i// l/ (8/xm)/^/ l/ (2.7/^m) > 3.16 at z ~ 2, and the 
mid-IR slope a > 2.1 if f v cx v~ a . Figure [T] il- 
lustrates how this color selection works with two 
simple example spectra. At z ~ 1 and z ps 2, the 
#(24, 8) ratio of a starburst ULIRG w ith strong 
PAH (NGC6240) (|Armus et al l 120061) is redder 
than that of an AGN (Mrkl014) \ Armus et"aT 



x For details of the XFLS observation plan and the data re- 
lease, see http://ssc.spitzer.caltech.edu/fls 
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2004). This is due to the fact that strong, broad 
PAH emission near 8 and 12 jxm moves into the 
24 (im filter at z ~ 2— 1. Obviously, this selection 
also picks up AGNs without much PAH emission 
but with steep red continua. 
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Fig. 1. — This plot illustrates how 24-to-8/im 
color can be used to select sources with PAH 
emission. It shows the mid-IR SEDs in the 
rest-frame wavelength in micron versus vj v for 
two ULIRGs which have been observe d by the 
Spitzer IRS (|Armus et all |2004 l2006j) . Here 
we extended the SEDs down to 2 micron us- 
ing the 2MASS an d HST/NICMOS (N GC6240) 
broad photometry ( Scoville et al. 2000). The Y- 
axis uf u is arbitarily scaled. The shaded re- 
gions show the widths of the 8 /im and 24 /im fil- 
ters at the rest-frame for z — 2. The shaded 
areas mark the rest-frame wavelength of the 8 
and 24 /zm filters at z = 2. The monochro- 
matic luminosity ratio between 24 and 8 micron 
(i.e~R(24,8) = vf v {2tym)/vf v (&iim)) is redder 
for starburst dominated ULIRGs like NGC6240 
than AGN dominated ULIRGs like MrklQ14 at 
z ~ 2. 

Photometry for all targets in the three filters, 
24 /im, 8 yum and R, are tabulated in Table [TJ The 
final photometry is slightly different from what 
was initially available for the GO-1 target selec- 



tion during the early period of the XFLS data re- 
duction (see Figure [9] in Section 4) . The change 
in MIPS 24 /im fluxes is very minor, 5% or less. 
The bigger change is in IRAC 8 /im fluxes, partly 
because of the new and more accurate aperture 
corrections, and partly because most of our tar- 
gets are very faint or not detected. This explains 
why i?(24, 8) values for some of the targets are 
slightly smaller than 0.5, a cut initially set by the 
early IRAC photometry. The full analyses of the 
selection function using Monte Carlo simulations 
will be performed in a separate paper (Sajina et 
al. 2006b). The brightness distributions at 24 /xm 
and i?-band for our target sample are shown in 
Figure O The majority of our targets are lmJy 
sources and the sample has a fairly narrow dy- 
namic range in 24 /im fluxes. The sources have 
R > 22, and many of them are fainter than the 
i?-band limit of 25.5 magnitude (shown as the last 
bin in the R distribution plot). The faint opti- 
cal magnitudes motivated in part this survey to 
measure redshifts using solely the mid-IR spectra. 
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Fig. 2. — These two panels show the distributions 
of 24 /im flux density and i?-magnitude for our 
sample. The last bin with the arrow in the R 
distribution plot shows the sources not detected 
in the R images, and fainter than the R limiting 
magnitude of 25.5. 
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2.2. The IRS Observations 

All observations were taken using IRS Staring 
mode, with nearby bright stars as peak-up targets 
to accurately center the slit on the science tar- 
gets. For the 52 targets, we obtain low resolution 
(^A- — 64 — 128) spectra in the long-low mod- 
ule of the IRS (LL, 14 — 40/im). For nine sources 
with IRAC 8 fim fluxes greater than 150/iJy, we 
also obtained Short-low, 1st order spectra (SL, 
7.5— 14/iirt). For IRS low resolution staring mode, 
one cycle consists of two exposures at two different 
positions on the slit (two nods). For our sample, 
one cycle of 242 seconds at each nod position is 
taken for the SL 1st order observations, and three 
to nine cycles of 122 seconds were used for the LL, 
depending on the brightness of the source. TableQ] 
shows the MIPS 24 /im ID, source ID named after 
the IAU convention, the broad band 24 /im, 8 /im 
and R-band fluxes for the 52 sources. The source 
ID following the IAU convention provides the in- 
formation on the source equatorial positions. We 
also listed the integration time of each individual 
exposure, the number of repeats at each nod po- 
sition for SL 1st order and LL 1st and 2nd order, 
which can be used to compute the total integration 
time for each order. For completeness, we include 
the initial ei ght sources whose spectra have been 
published in lYan et al ] (l2005r ). Our targets are 
all unresolved in all IRS wavelength. Since the 
IRS slit width is 3.7" and 10.7" for the SL and 
LL respectively, and the full-width-half-maximum 
(FWHM) at 8 /mi and 24 /im is ~ 2" and ~ 6", 
light loss correction due to the finite IRS slit width 
is not needed for our observations. 

3. Data Reduction and Analyses 

3.1. Two-Dimensional (2D) Spectral Imag- 
ing Processing 

Our reduction starts with the 2D Basic Cali- 
brated Data (BCD) produced by the IRS pipeline 
(S13 version) at the Spitzer Science Center (SSC). 
The processing steps taken by the IRS pipeline 
include ramp fitting, dark sky subtraction, droop 
correction, linearity correction, flat fielding, and 
wavelength and flux calibration. We perform ad- 
ditional processing, including subtraction of the 
"residual" background in the BCD images, mask- 
ing any significantly deviant pixels, and stack- 
ing images taken at the same nod position. We 



build the background image by using only the 
BCD images of the same target, including both 
nod positions. At the LL, we have roughly (12- 
32) BCD images for computing the background. 
Our method is iterative. In the first pass, a me- 
dian combined background image is obtained after 
masking out the target spectra. After subtracting 
this background image from each BCD image, we 
go through the second iteration of identifying and 
masking out serendipitous spectra in the images, 
then produce the final background estimate. To 
minimize the effects of the deviant pixels, we used 
a biweight estimator H to compute the background 
and its dispersion. The noise in each pixel of the 
final images is computed by adding quadratically 
the noise evaluated from the background to the 
Poisson noise from the target spectra. The er- 
ror spectrum is measured from the corresponding 
noise image. 

Before spectral extraction, we identify and cor- 
rect the "rogue pixels" in the 2D images. The 
"rogue pixels" are pixels whose dark current is 
abnormally high and varies with time and differ- 
ent sky background radiation. These pixels are 
highly deviant from the local background. We 
identify the "rogue pixels" by computing the dis- 
persion of the noise around every pixel and flag- 
ging pixels which are above 5er of the mean value. 
This threshold is relatively conservative. For any 
deviant pixel which does not fall on the spectra 
(within one FWHM profile), the code automati- 
cally sets its value to zero (which is the value after 
background subtraction). For the deviant pixels 
which do fall on the target spectra (within one 
FWHM), we visually examine each pixel, and re- 
place its value with the linear interpolation (±1 — 2 
pixels) along the wavelength direction. The inter- 
polation is appropriate for our data because most 
of the detected spectral features (i.e. PAH) are 
much broader than 1-2 pixels, and will not be 
smeared significantly. This type of correction is 
usually needed for only one or two pixels per spec- 
trum. 



Reference: "Data Analysis and Regression: A Second 
Course in Statistics" , Mosteller and Tukey, Addison- 
Wesley, 1977, pp. 203-209. 
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Fig. 3. — The mid-IR spectra for 47 sources with redshift measurements. The spectra are presented in 
observed flux density f v in mJy versus the rest-frame wavelength in micron. The spectra were smoothed by 
a 3 pixel boxcar, in order to enhance the broad features such as PAH emission and silicate absorption. In 
each panel, the redshift, source ID and the spectral type (see text for the description) are labelled at the 
top-left corner. The large, open squares indicate the broad band, 24 /im photometry, and the gray shaded 
region marks the ±lcr error for the spectrum. The error spectrum is not smoothed. The error computation 
is discussed in §3.2. The dashed lines mark the PAH emission features, and the arrows mark the silicate 
absorption. 
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Fig. 3. — Continue. 
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Fig. 3. — Continue. 
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Fig. 3. — Continue. 
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Fig. 3. — Continue. 
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Fig. 3. — Continue. 
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Fig. 4. — The five observed spectra that did not yield any redshift measurements. The spectra are plotted 
in the observed flux density versus observed wavelength. As in Fig [3l the open squares indicate 24 //m broad 
band photometry. The five spectra have power-law continua without identifiable spectral features. The 
spectral type is set to t = 4. 
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3.2. Spectral Extraction 

Using the steps described in the section above, 
a final coadded 2D spectral image is produced 
for each nod position for each IRS slit. One- 
dimensional (ID) spectra were extracted, then av- 
eraged to produce the final ID spectrum for each 
module. The spectral extraction is done using the 
SSC software SPICE vl.Sbeta version^. In SPICE, 
the spectral extraction aperture can be adjusted to 
maximize the S/N ratio. For our data, we used 
an extraction aperture of 6 pixels at 27 /im 1 pixel 
narrower than the default PSF size in SPICE, and 
the shape of the extraction aperture simply fol- 
lows the PSF size as a function of wavelength. 
This method minimizes the uncertainties due to 
the background, and maximizes the S/N ratios for 
faint spectra. We used the simple method of sum- 
ming up the light within the extraction aperture. 
The flux scales of the final, averaged SL and LL 
spectra are consistent with the IRAC 8 /jm and 
MIPS 24 /im broad band fluxes. No rescaling is 
applied to any of the spectra. The final spectrum 
for each source is produced by averaging the spec- 
tra from the two nod positions and combining the 
spectra from the 1st, 2nd and 3rd order together 
using linear interpolation. Different orders have 
different pixel-to-wavelength scales, therefore, we 
resample the overlapping spectra onto the small- 
est scale, and produce the final combined spectra 
by averaging. We found that the flux calibration 
between different orders is consistent. 

4. Results 

4.1. Redshift Measurements from Mid-IR 
Spectra 

The observed 52 spectra from the sample are 
presented in Figure [3] and Figure [4] Of the 52 
sources, 47 have redshifts measurable from their 
mid-IR spectral features, such as PAH emission 
and silicate absorption, yielding redshift efficiency 
of 90%. The remaining five sources have power- 
law spectra without any identifiable features, thus 
have no redshift measurements (Figure |4|). For 
the majority of our spectra, the redshifts are de- 
termined by the multiple spectral features, and 
are thus secure. As in any spectroscopic redshift 



3 see http://ssc.spitzer.caltech.edu/postbcd/spice.html for 
the description of the software and the public release status. 



studies, for a few sources, the spectral identifica- 
tions could have some ambiguities. We therefore 
identify our redshift measurements by a quality 
flag Q z , with Q z = a for secure identification, and 
Q z — b for redshifts spectral features whose iden- 
tifications are ambiguous. We have only 6 z with 
Q z = b. If excluding these 6 redshifts, our red- 
shift efficiency is 80%. Flag Q z and redshift errors 
are two different parameters which quantify the 
security of the spectral identification and redshift 
accuracy respectively. 

The redshift errors are estimated as follows. We 
compute individual (l + z); = X bs(i)/ X re st{i), us- 
ing spectral feature i, based on the measured fea- 
ture centroid X bs(i) and the corresponding rest- 
frame wavelength X rest (i). The final redshift z is 
the average value of Zi, and the la deviation is 

i z . ~ < ^ z ~ > \2 

computed as la — — % J n _ x ■ Another source 

of redshift errors is Azi, from the centroid un- 
certainty AXi for each feature. The total z error 
should be the quadratic sum of the la standard de- 
viation and Azi . For spectra whose redshifts based 
on a single spectral feature, the uncertainty is only 
estimated from the uncertainty of the feature cen- 
troid. For example, for spectra with only silicate 
absorption, the redshift errors can be as large as 
0.1-0.2 because of the intrinsic broad width of 
this feature. For spectra with both PAH emission 
and silicate absorption features, we use only PAH 
features for measuring redshifts since they are gen- 
erally narrower than silicate absorption, thus give 
more accurate results and smaller redshift errors. 

One of the critical parameters in determin- 
ing the redshifts are the rest-frame wavelength 
Xrestii) for various PAH features. These broad, 
complex emission bands roughly have similar rest- 
frame centers, although small variations have been 
found in studies of local PDRs (Vermeij et al. 
2002; Peeters et al. 2004). For our purpose, 
we measure X re st(i) from the mid-IR spectrum 
of a local starburst galaxy NGC7714, observed 
by IRS in the same config uration and resolu - 
tion as used for our sample (jBrandl et al.l l 20041) . 
Taking into account the distance, we computed 
the rest-frame wavelengths of PAH emission fea- 
tures based on the measured centroids. We used 
the following values for the reference wavelength 
X rest (i), 6.2180 ± 0.0061, 7.7118 ± 0.0229, 8.5972 ± 
0.0077, 11.2665± 0.00353, 12.8030 ± 0.0077 ^m for 
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6.2,7.7,8.6,11.2 and 12.7/zm PAH features respec- 
tively. 

Table [2] lists the redshift, the redshift error, 
the estimated rest-frame luminosity at 5.8 fim, 



vL v (5.Bix) 



the redshift quality flag 



(Qz), the spectral group classification (t), and the 
spectral features used for redshift measurements. 

We obtained optical and near-IR spectra for a 
subset of the IRS sample using LRIS, DEIMOS 



and NIRSPEC on the Keck telescopes (|Oke et al 



1994 iFaber et al.ll2003t lMcLeanll2000l ). We were 
able to determine 17 redshifts from the Keck spec- 
tra. All of the Keck redshifts confirm the estimates 
based only on the mid-IR features, even in a few 
ambiguous cases where the IRS redshifts are based 
only on weak features. Figure [5] illustrates the 
Keck z versus the IRS-based z. For two sources 
with weak features, MIPS110 and MIPS133 with 
t = 3 and Q z = b, their IRS redshifts are primar- 
ily based on the mid-IR continuum slope change 
around 10 /im. These results have been confirmed 
by their Keck spectra, suggesting that our visual 
identification of weak mid-IR features is reason- 
able. This 10 /im break is probably not a new fea- 
ture, but due to small amount of silicate absorp- 
tion at 9.8 /im, and generally the slightly curved 
shape at the break supports this idea (see Fig- 
ure [3|). One of the five sources without IRS- 
based redshifts, MIPS279, was observed with the 
Keck LRIS, and the redshift is at 1.23. Tabled 
shows the redshifts from the Keck and the IRS. 
We also listed the source type classified from the 
IRS, as well as from the optical spectra when it 
is possible. At the near-IR, we have only enough 
wavelength coverage to detect H a , thus definitive 
source type classification is not possible. However, 
all H a emission lines are narrower than 800 km/s. 
As shown, the redshift and the classification from 
the optical/near-IR data are consistent with what 
measured from the Spitzer spectra. 

4.2. Mid-IR Spectral Type 

Examining Figure [3] and Figure |4j it is ap- 
parent that the 52 spectra in our sample fall 
broadly into three categories: strong PAH emis- 
sion sources, sources without PAH but with deep 
silicate absorption at 9.8 /im and mid-IR contin- 
uum dominated sources with weak or no emis- 
sion/absorption features. Making this crude clas- 
sification according to broad spectral features al- 



lows us to have a quick assessment of the main 
properties of the sample. PAPER II of this se- 
ries discusses in detail the full spectral fitting 
and analyses. We use numeric t = 1 to designate 
systems with strong PAH emission. Quantita- 
tively, the spectra classified by eye in this cat- 
egory roughly correspond to 7.7 /im rest-frame 
equivalent width (EW) greater than 0.7 /im, or 
equivalently, 6.2 /xm EW„ S ( > 0.2 /im (computed 
in PAPER II). Sources without any evidence for 
PAH emission, but strong 9.8 /mi silicate absorp- 
tion feature are designated t — 2. Intermediate 
sources between these two types, i.e. with strong 
Si absorption, as well as prominent 8 /im bumps, 
which are possibly due to PAH emission (occa- 
sionally supported by 6.2 /im features) are desig- 
nated t = 1.5. The rationale for separating this 
sub-group (t = 1.5) is that a good fraction of high- 
z, IR-luminous sources seem to have this type of 
spectra, wea k PAH emiss i on wi t h strong silicate 



absorp tion ( Houck et al. 2005; IWeedman et al 



2006aJ). The mid-IR continuum dominated galax- 



ies are separated into sources with strong con- 
tinuum plus weak PAH emission and/or silicate 
absorption features (t = 3), and sources with pure 
power-law continua without identifiable features 
(t = 4, thus without redshifts). As shown in Fig- 
ure [3l we have 17 sources with strong PAH emis- 
sion (t = 1), 17 sources with only deep silicate ab- 
sorption (12 with t = 1.5, 4 with t — 2), 13 spectra 
with i = 3 whose emission and absorption features 
are weak, and finally, 5 sources with power-law 
continua without redshifts. Some of these five 
sources could be AGNs with very red mid-IR con- 
tinua, and the rest could also be starbursts at 
z > 3, where most strong PAH emission and sili- 
cate absorption features have redshifted out of the 
IRS longest wavelength window (5 — 38 /im) . 

4.3. Redshift distribution versus spectral 
type 

Figure [6] presents the redshift distribution of 
all 47 redshifts, as well as the z-distribution for 
each source type. Of the 47 redshifts, 12 are 
0.6 < z < 1.3 and the remaining 35 are beyond 
z = 1.5. This indicates that our selection cri- 
teria are indeed very effective in selecting z > 
IR luminous galaxies. The \ow-z contamination 
comes from the fact that at 0.6 < z < 1, IRAC 
8 /im and MIPS 24 /im filters sample the rest-frame 
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Fig. 5. — We compare the redshifts measured from 
the Keck optical and near-IR spectra with the IRS- 
based redshifts. The straight line marks the 1 : 1 
ratio. 



4 /tm continuum and 12.7/tm PAH respectively. 
The strong 12.7 (im PAH emission in these low- 
z sources boosts the 24/8 /tm color into our red 
color cut. The dip at 1.3 < z < 1.6 in Figure [5] 
is due to the fact that at z — 1.5, the 24 /mi filter 
samples the rest-frame 9.8 /im, thus any samples 
with 24 /im flux cutoff will inherently bias against 
selecting sources with 9.8 [im silicate absorption 
around redshifts of 1.5. Sources without strong sil- 
icate absorption should be selected. For example, 
we have two sources, MIPS15928 and MIPS8342, 
at z = 1.52 and 1.57 in our sample, which have 
weak or no 9.8 /tm absorption as shown in Figure[3] 
As shown in Figure all of the 16 sources with 
only deep silicate absorption (t = 1.5 and t = 2) are 
at redshifts beyond 1.5. Of the 17 sources with 
strong PAH (t = l), 7 are at z ~ 0.61 -1.0, and 
the remaining 10 are at z ~ 1.52- 2.59. For the 
13 sources with t = 3, 6 are at z <~ 0.7- 1.23, and 8 
at z~ 1.57-2.4. 

4.4. Stacking the mid-IR Spectra 

To increase the S/N ratio and to examine the 
variation and uniformity of the spectral features, 



Fig. 6. — Redshift distributions separated by the 
spectral types. The black histogram is for all of 
the redshifts (47), and the red is for 17 of the 
strong PAH sources (t = l), the green for 17 of 
the highly-obscured sources (i = 1.5& 2), the blue 
for 13 of the continuum dominated objects with 
weak features (emission and/or absorption). 

we applied stacking analyses to the spectra sep- 
arated by each class and by redshift (z < 1 and 
z > 1). The spectra are normalized either at 4.8- 
5.2 /tm for sources at z > 1, or at 15 /im for z < 1. 
We computed the median spectra as well as 20 and 
80 percentile spectra to illustrate the spectral vari- 
ation. For the strong PAH class (t = 1), because of 
different rest-frame wavelength coverage, the spec- 
tra at z ^ 1 and at z ^ 1 are stacked separately. 
The same approach is taken for type t — 3, except 
that we exclude two spectra in the z > 1 bin due to 
their much steeper continua. The stacking results 
are shown in Figure [71 The improved S/N median 
spectra confirm the spectral features identified in 
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individual spectrum. Particularly, for t = 3, z > 1 
spectra, although the PAH features in individual 
spectrum are quite weak, their presence is indeed 
confirmed by the median stacked spectrum for this 
class. The median spectrum for t = 3, z < 1 ob- 
jects shows the weak absorption trough at 9.8 /im, 
although truncated, confirming the redshift deter- 
minations for this group of sources. 

Similar types of mi d-IR spectra have b een seen 
among local ULIRGs (jArmus et al.ll2006l ). In Fig- 
ure [3 we overlay in red lines the rescaled, IRS 
spectra of local ULIRGs. We emphasize that these 
local spectra were plotted not for the purpose 
of spectral fitting, but to illustrate their similar 
shapes to our spectra at z ~ 1 — 2. For exam- 
ple, our type t = 2 spectra are similar to the local 
ULIRG IRAS F 08572+3915, which is a highly em- 
bedded system (jSpoon et al.ll2005| ). Type t = 1.5 
spectrum has a broad peak around 7 — 8 /tm. The 
closest example we can find among th e local, avail- 



able ULIRG spectrum is Arp 220 (jArmus et al 
2006). It is possible that some of that broad 
7 — 8 /mi peak could be from PAH emission, and 
not just simple mid-IR continuum. For t = 3 spec- 
tra, we plot the pure power-law continua, indi- 
cated with red lines in the figure. The slopes (a, 
f u oc v~ a ) are 2.2 and 1.8 for the median spectra 
at z > 1 and z < 1 respectively. These slopes are 
much steeper than the mid-IR spectra of classic 
AGN, like NGC1068, with a~l. 

Figure [7] illustrate three broad categories of 
spectra observed in our sample, spectra with clear 
PAH emission, ones with only deep silicate ab- 
sorption and continuum dominated spectra with 
only weak PAH and/or silicate absorption. Tak- 
ing the median spectra, we measured the optical 
depth at 9.8 /tm, r 9 . 8Atm = ln(S , cont /<S'min), by sim- 
ply approximating the continuum with a power- 
law between 7 /tm and 13 /tm. This crude esti- 
mates show silicate absorption is significant for 
both strong PAH sources (t = 1) as well as absorp- 
tion only systems (t — 1.5,2). These spectra are 



optically thick at 9.8 /jm with Tg. 



1.1 - 1.3, 



implying that close to 60% of our sample are heav- 
ily obscured at mid-IR. We note that for systems 
with strong PAH emission, the measured optical 
depth at 9.8 /im silicate absorption provides only 
a lower limit since the broad wings from 7.7, 8.6 
and 11.3 /im PAH emission tend to quickly fill in 
the absorption trough, resulting under-cstimate of 



the true opacity. 
4.5. Luminosities 

4-5.1. 5.8 [im monochromatic continuum lumi- 
nosity 

To characterize the energetic output of our 
sources, we measure the rest-frame 5.8 /tm monochro- 
matic, continuum luminosity, Ls.g = pL u (5.8^,m) 
from our spectra. This is the rest-frame wave- 
length where we have data for most of our sources. 
Table [2] includes the derived logLs.s luminosities 
for our sample, which range between ~ 10.3 for the 
lowest-z sources up to ~ 12.6. Figure [5] shows the 
distribution of L5.8 versus z for our sample with 
sources with and without strong PAH emission 
in red and blue symbols respectively. We found 
that some of the strong PAH emitting galaxies 
at z ~ 2 are ~0.5de x more luminous th an lo- 
cal starburst ULIRGs (jArmus et al.ll2006l) . 



The 

mid-IR luminosity trend with redshift is essen- 
tially our selection function. The minimum and 
maximum 24 /im fluxes are translated to L5.8 us- 
ing NGC6240 spectrum, shown as the dashed 
lines in the figure. The peak and the valley at 
z ~ 1.5 and z ~ 2 reflect the silicate absorp- 
tion trough and the strong 7.7 /tm PAH emission 
passing through the observed 24 /tm filter at these 
redshifts. Studies of local ULIRGs suggest that at 
Lir, = Ls-iooo^m ~ 10 12,5 L(T) ULIRGs tend to be 
AGN -dominated ( Lutz et al. 1998; Vcill eux et al.l 
Il999l ) . In Figure we do not observe any sig- 
nificant trend in mid-IR luminosity versus source 
type. This does not imply that our high-z sam- 
ple is different from the local ULIRGs. To study 
the luminosity versus source type, we need to de- 
rive acc urate Ltr using long wavelength data at 
1.2mm ([Lutz et al.l l2005h and MIPS 70, 160 /tm. 
We defer this analysis to a later paper. 

However, in order to put our sources in con- 
text with other galaxy populations, here we use 
a crude estimate of Lm based on £5.8- This 
is necessarily very crude, as 5.8 /tm does not di- 
rectly sample the far-IR thermal emission domi- 
nating the bolometric infrared emission. Brandl 
et al. (2006) showed that using the 15 /tm and 
30 fim fluxes together addresses both these issues, 
and provides a good conversion to Lir. For the 
z<1.6 sources (about 1/3 of the sample), using 
both the IRS spectra and the MIPS70/tm points 
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Fig. 7. — The six panels snow the median spectrum (in black lme) per type for z ^ 1 and z > 1 sources. For 
type t = 1 and t = 3, we stacked spectra for z < 1 and z > 1 separately. For class t = 3 and z > 1 bin, we 
excluded two spectra in the stacking (see the text for detail). Type t — 1.5 and t — 2 have only z > 1. The 
cyan shaded regions indicating the 20 and 80 percentfile. The source type, average z, and number of spectra 
included in the stacking are indicated at the top-left corner of each panel. The red spectra in the first four 
panels represent the scaled local ULIRG spectra to illustrate the similarity between spectra at high-z and 
2 = 0. The local ULIRG spectra were chosen from Armus et al. (2006b) to have the closest match with 
the median spectra. In the last two panels, for t — 3 class, the red lines indicate power-law continua with a 
(f v oc v~ a ) of 2.2 and 1.8 for z > 1 and z < 1 respectively. The slope a is much steeper than that of the 
local AGN spectra like NGC1068. 
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we can derive both these points (see PAPERII 
for details). Using these and the Brandl et al. 



relation, we find and average Lir/Ls., 



25 



+10 

-8 



This average ratio is consistent with the value de- 
rived from the IRS spectra of the local ULIRGs 
sample, except for the extremely cold and warm 
outliers Arp220 and IRAS08572+3 915 (for this 



sourc e, the ratio is the lowest, 5) (jArmus et al 
2006). Given the 5.8 /im luminosities in Table 2 



and the above ratio implies that our sources have 
L m ~ 10 12 - 10 13 L Q and are ULIRGs or type II 
QSOs. With the average of < L 5 . s > = 9xlO n L Q , 
even with the lowest ratio of 5, the implied Lib is 
still 4.5 x 10 12 i© = 2 x 10 46 ergs/s, similar to the 
energy output from quasars. 
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Fig. 8. — Monochromatic, rest-frame 5.8 /im con- 
tinuum luminosity vs. redshift with different sym- 
bols marking the different source types. The red 
and blue color crudely separates the sources with 
strong and weak PAH emission respectively. The 
dashed line shows the L 5 s vs. z relation if we take 
a SED like NGC6240 with the F obs (2A^m) scaled 
to the lower and upper limits of our survey sam- 
ple. The errorbar shows the mean errors for our 
sample. 

4-5.2. The PAH fraction of mid- 1 R luminosity 

In normal star- forming galaxies, PAH emission 
accounts for ~ 10-15% of th e total infrared lumi- 
nosities ([Smith et all [2006). For local ULIRGs, 



crude estimate by measuring the PAH fraction 
in the median spectra with strong PAH (see 
Fig. [7]). A more careful spectral fitting of in- 
dividual spectrum is done in PAPERII. We es- 
timate the fraction of summed PAH luminosity 
over 5.8 /im, monochromatic continuum luminos- 
ity, LpAn/i/L v (5.8fj,m). The sum of the PAH 
luminosities, ^Lpah, includes the 6.2, 7.7, 8.6 
and 11.3 /im features. We use a simple straight- 
line continuum between 5 /im and 14 /im, which 
allows for a Si absorption feature at 10 /im. The 
PAH flux is estimated by integrating the full 
features from the continuum-subtracted spec- 
tra. Given the 20th and 80th percentile range 
in the median spectra shown in Fig. we ob- 
tain X) l pah/^ M5.8/im) ~ 0.2-0.7. This im- 
plies that the PAH emission accounts for ~ 1 - 3% 
of the total infrared luminosities (£ ir), which is 
comp arable to that of local ULIRGs (|Armus et al. 
2006)) , but lower than pure starburst galaxies from 
Brandl et al. (2006). 

This is much smaller than what is measured 
among normal galaxies, and similar to that of local 
ULIRGs. 

4-5.3. The rest-frame UV luminosities 

Lastly, we examine the rest-frame UV versus 
IR luminosities, ii60o/^iRj which indicates how 
much UV optical photons being absorbed and re- 
processed by dust, thus the averaged dust obscu- 
ration. We compare our sample with other well- 
known UV-selected or IR-selected high-z popula- 
tions. The derivation of log Lib for ° ur sample 
was discussed in § 14.5.11 We derive the UV lu- 
minosities from the observed i?-band photometry 
and redshifts, by assuming a flat UV spectrum 
(i.e. uF„ = const, /3 = -1, for A<4000A). This 
assumption of /? = — 1 is not too far off since SMG 
and LBG ha ye the median ' value of -1.5 and -1 
respe ctively ([Chapman et al. I l2004bt iReddv et alj 
2006). In addition, since at z ~ 2, the fi-band 



this fraction is small, on the order of a few per- 
cent (Armus et al. 2006b). Here we give a 



samples the rest-frame 2333A, the K-correction 
in the estimate of £i600 from i?-band photome- 
try should be small. Most of our sources have 
£ieoo £3x lO 10 ^©. Figure [TO] shows L m /L 160Q 
vs. Lm + ^1600 f° r our sample, in comparison with 
other galaxy samples compiled by Reddy et al. 
(2006). The lower luminosity, lower redshift sub- 
sample of our galaxies overlap the most-luminous 
LBGs (a.k.a. UGRs) and the BzKs. The bulk of 
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our sample however are closer to SCUBA galax- 
ies, although our galaxies are often more luminous. 
The mean Lih/Liqqo ratio is about 5 for local nor- 
mal galaxies, and much higher for high-z galax- 
ies, with (1-100) for LBGs+BzKs and 100 - 3000 
for SMGs and our sample. Figure [TO] suggests 
that more luminous galaxies tend to have higher 
IR/UV luminosity ratio, and on the whole our 
sources are among the most luminous and most 
obscured sources known. In particular, the local 
galaxy sample (open squares, Bell 2003) and high- 
z galaxies seem to follow different correlation be- 
tween L to t vs. Liji/Liqoo (dashed and solid lines). 
At a fixed IR-to-UV luminosity ratio, high-z galax- 
ies are more luminous than z = counterparts, 
and at a fixed total luminosity, high-z galaxies 
have lower -Zj/_r/Li600 ratio than that of local ones 
(although local galaxies with very high IR-to-UV 
ratios are very few). 

5. Discussions and Summary 

5.1. Nature of the PAH emitting ULIRGs 

at z ~ 2 

The detection of PAH molecules in one third 
of the sample provides the direct evidence that 
aromatic features are already an important com- 
ponent of galaxy spectra at the early cosmolog- 
ical epoch. Quantitatively, the total energy in 
PAH features is roughly 20-70% of the 5.8 /im 
continuum luminosity (^L„(5.8/im)), correspond- 
ing to ~ 1 - 3% of the total infrared luminosity. 
Studies utilizing IR SEDs need to properly con- 
sider the strength and distribution of PAH fea- 
tures since they have significant impact on the 
model predictions, such as redshift distribution 
and l uminosity functions of 24 jxm selected sam- 



ples ([Lagache et al 



20041 ) . In addition, the exis- 
tence of PAH molecules in z ~ 2 ULIRGs implies 
that metals from massive young stars formed at 
earlier epoch have gone through a complex pro- 
cess — cycling through dense molecular gas, form- 
ing photo-dessociation region (PDR) and starting 
another intense star formation again. 

The estimated rest-frame 5.8 /jm luminosities 
for our strong PAH sources are a factor of 5- 
10 times more luminous than the local starburst 
dominated ULIRGs. There are two possible phys- 
ical explanations for this luminosity increase. The 
first scenario is that these z ~ 2, strong PAH 



sources have much higher star formation efficiency 
as well as higher gas-to-mass fraction than lo- 
cal ULIRGs. The efficient star formation for- 
mation converts the large gas reservoir to stars 
and produce higher luminosity. This hypothesis 
is supported by studies of z ~ 2 sub-mm galax- 
ies, which have found that SMGs are scaled-up 
version of local ULIRGs, in the sense of lumi- 
nosity , mass, star formation rate and ga s frac- 



tion (iChapman et all l2004ri iGreve et all 120051: 



Tacconi et al.l 2006). It is possible that the z ~ 2 
Spitzer ULIRGs and SMGs are like local ULIRGs, 
and are in the stage of "maximum starburst" , 
i.e. an extreme mode of star formation where gas 
mass density is so high that it is balancing out the 
momentum-driven wind (radiative pr essure, su- 
pernova explosion and st ellar winds) (jElmegreen 
19991 : iMurrav et al.l liooil ). resulting in the most 
efficient star formation. The second possible ex- 
planation is that although these z ~ 2, PAH emit- 
ting ULIRGs are dominated by starburst at the 
mid-IR, they could still have higher AGN contri- 
butions than the local starburst ULIRGs. Our 
spectra, with low spectral resolution and limited 
wavelength coverage, can not revealed the AGN 
observational signatures, such as high ionization 
ionic lines. 

5.2. Nature of the non-PAH sources 

The AGN unification model as well as the ob- 



served hard X-ray background (jAntonuccil 11993) 
have predicted the existence of type II quasars, 
with a rough ratio of type II to I of (4-10):l 



(|Madau et al.l 11994k IComastri et al.lll995h . How 



ever, UV and X-ray obseryation s have so far failed 
to find many (jStern et alj|2002t) . Our AGN dom- 
inated sources have L 5 , s > 9 x 1O 11 L = 3 x 
10 45 ergss -1 , corresponding to a total infrared lu- 
minosity on the order of 10 47 ergs s _1 . Such a high 
luminosity is comparable to the bolometric lumi- 
nosity of bright quasars. If most of this energy is 
generated by dust re-processing of the UV, soft X- 
ray photons from a black hole accretion disk, the 
inferred black hole mass is on the order of 10 8 M Q , 
assuming Eddington accretion rate. Even though 
one third of the sample do not show much mid- 
IR absorption along the line-of-sight, all of our 
sources are surrounded by substantial dusty ma- 
terial, as indicated by their large rest-frame IR to 
UV luminosity ratio. The deep silicate absorp- 
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tion systems probably have dust grain distribu- 
tion with high filling factor, whereas the objects 
with weak absorption features have lumpy dust 
distributions. Detailed modeling of deep silicate 
absorption spectra by Levenson et al. (2006) sug- 
gests that these sources contain central nuclear 
sources deeply embedded in a smooth distribution 
of dust that is geometrically and optically thick. 

Enormous infrared luminosity and high dust 
content suggest that our AGN sources are poten- 
tially type II quasars or at earlier stages before be- 
coming supermassive black holes. Of the total 52 
sources, 22 are AGN dominated systems at 1.75 < 
z < 2.7. This gives the co-moving volume density 
of 5x 10~ 7 Mpc~ 3 , for the survey area of 3.7 sq. de- 
gree. The type I quasar co-moving space density at 
z ~ 2 is about 6 x 10~ 7 ' Mpc~ 3 with vL v {%jim) > 
4.1 x 10 45 ergss" 1 - 10 12 L Q (|Brown et al.ll2006l ). 
This luminosity limit is roughly comparable to the 
averaged L5.8 of our sample. The space density 
measured from our sample is only a lower limit to 
all of type II quasars. At face value, these space 
density estimates are consistent with the specula- 
tion that our AGN dominated systems could be in 
the early stages of quasars evolution. 

The AGN sources in our sample also have very 
red mid-IR (5-15 ^m) slope, with a > 2 (f v oc 
v~ a ), much redder than the averaged slope of 1.3 
for Palomar-Gree n (PG) quasars at z ~ 0—1.2 
(jHaas et al.ll2003l ). The IRS spectra of our sample 
demonstrate that mid-IR colors (rest-frame, broad 
band) can not separate high-z PAH emitting star- 
bursts from AGNs with red mid-IR continua. For 
example, studies using blue 70 /im-to-24 (im color 
to select AGNs will miss a population of AGNs 
with red mid-IR continua, similar to our sources. 

5.3. Relation with other high-z galaxy 
populations 

We compare our sample directly with the Houck 
et al. (2005) Spitzer selected z ~ 2 sample. 
The IRS GTO team obtained low resolution IRS 
spectra for 58 sources with / 1/ (24/xm) > 0.75 mJy 
and / > 24mag (vega) in the Bootes field. Al- 
most all of their observed spectra have deep sili- 
cate absorption/pure continuum only and no ob- 
vious, strong PAH emission, e xcept a few spectra 
with broad peaks at 7 — 8^m ( Houck et al.ll2005t 
IWeedman et al.ll2006bh . Of their 58 sources, 76% 
(44) have redshifts based on the IRS spectral fea- 



tures and the remaining 24% are pure power-law 
without yielding any z measu rements. The re- 
sults from the IRS GTO study (|Houck et al.ll2005t 



Weedman et al. 2006b) differ from ours in the 



sense that they detect much less fraction of sources 
with strong PAH emission. This difference can be 
explained by the difference in the target selections. 
The IRS GTO study uses only two criteria: 1). 
iWm > 0.75mJy, 2). v ] v {2^m) j v j v (l) > 100, 
and a few sources have v f v {2^m) / v f v {I) > 75 — 
100. 

Figure [9l in three panels, to illustrate 24/8 
and 24/R color-color diagram and the number of 
sources as functions of the two colors. Here the red 
and blue colors indicate objects with and without 
strong PAH respectively. The color is computed 
as vf v ratio between two wavelengths. Both Panel 
(a) and (b) show that sources without strong PAH 
tend to have redder 24/R colors (blue lines) in 
comparison with sources with strong PAH (red) . 

The Houck et al. sample has much redder (more 
than a factor of 10 in flux ratio) mid-IR to opti- 
cal colors than our sample, occupying the regions 
redder than the shaded area in Panel (a), which 
indicates the GTO 24/im-R flux ratio limit. Here 
we convert the GTO I magnitude to R used by our 
sample assuming vf v oc constant (flat spectrum). 
Most of our PAH strong sources lie the left side 
of the shaded area. The extremely high IR-to- 
optical ratio selects galaxies whose central energy 
engine is completely obscured by dust, very little 
UV photons escape from the centers, and further- 
more, there is very little star formation outside 
the nuclear regions. In contrast, some of our PAH 
emitting galaxies do have star formation outside 
of heavily obscured central regions. 

What about our sample in relation with other 
optical UV and (sub)-mm selected sample? The 
z ~ 2 ULIRGs revealed in our study represent 
an IR luminous population, which have very dif- 
ferent dust and energetic properties from those 
of SMGs and LBGs. Spitzer observations indi- 
cate that SMGs tend to have much fainter 24 /im 
fluxes (on average a few 100/iJy) and their mid- 



IR spectra have strong PAH emission (jLutz et al 
2OO5) ; IPope et al ' " 



2006). This is expected since 



Spitzer 24 fim selection tends to bias toward mid- 
IR bright, warmer, AGN type sources whereas 
the sub-mm selection prefers galaxies contain- 
ing substantial cold gas emitting at the far-IR 
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Fig. 9.— Panel (a) shows the 24/R and 24/8 /mi 
colors for the 52 sources we obtained the IRS 
low resolution spectra. The arrows indicate the 
sources with flux limits in R and/or 8/iim band. 
The red plus symbols indicate the sources with 
strong PAH (t = 1), and the blue open circles are 
for the deep absorption only spectra and mid-IR 
continuum dominated spectra (t = 1.5,2,3). The 
shaded region marks the 24/R color cu t used by 
the IRS GTO study (|Houck et al.ll2005l) . The ma- 
jority of the GTO sources have 24/R color redder 
than the color marked by the shaded area. Panel 
(b) and (c) are the histograms of the source distri- 
butions for 24/R and 24/8 /im respectively. The 
red and blue lines are for sources with and without 
strong PAH emission respectively. 



( Blain et all [2004). The initial 1.2mm obsesr- 



vations of a subset of our sources by Lutz et 
al. (2005) also suggests that Spitzer z ~ 2 
sources have on average less cold dust emission 
than SMGs. In addition, our z ~ 2 sources have 
very little overlap with UV-optical selected Ly- 
man break galaxies (LBG), which a re much less 



infrar ed luminous than our sources (jReddv et al 
20061 ). Of the sample of LBGs [z ~ 1.5 - 2.6) in 



the HDFN, 60% have 24 /im counterparts with the 
average flux de nsity of (20-30) u,Jy, corresponding 
to 5 x lO n L (|Reddv et al.ll2006h . Deep 24 /mi 



images reaching micron- Jy limits will probe the 
same galaxy population as SMGs and will have 
significant overlap with LBGs. 

We make crude estimates of space density of 
bright 24 /^m sources at z ~ 2, in comparison with 
other types of galaxies. Of all the 24 /im sources 
brighter than 0.9mJy over 3.7 deg 2 in the FLS, 5% 
(59) meet our sample selection criteria (our sam- 
ple has 52 of these 59) . Of the 52 sources from our 
sample, we have 34 at z ~ 1.5 - 2.7. Using this frac- 
tion to scale to the total 59 sources, the comoving 
space density of z ~ 2 ULIRGs satisfied our selec- 
tion criteria is n~ 2 x 10~ 6 Mpc~ 3 for z = 2.1±0.6 
and I/ir > 5 x 10 12 L Q (see Figure [8] for the limit). 
This volume density is only for ULIRGs with very 
red mid-IR color , and obviously is only a lower 
limit for the total population of z ~ 2 infrared 
luminous galaxies. For comparison, the comov- 
ing number density of sub-mm detected galax- 
ies at (z) = 2.2 is roughly 6 x 10~ 6 Mp c~ 3 for 



■Ltr > 4 x 10 12 L g ([Chapman et alJ l2004aj ). This 



implies that our survey has revealed a popula- 
tion of z ~ 2 ULIRGs which is as numerous as 
SMGs, and the total infrared luminous popula- 
tion at z ~ 2, including both ours (warm) and 
SMGs (cold), is at least 10 _5 Mpc~ 3 for L m > 4- 
5x 1O 12 L . 

The rest-frame UV color selected galaxies have 
co-moving space density of 2 x 10 _3 Mpc~ 3 for the 
B aimer br eak selected galaxies at (z) — 1.77 and 



(z)=2.32 ijAdelberger et al.ll2005h . Optical/UV 



selected galaxies are more numerous than IR lu- 
minous galaxies (Spitzer + SMGs) by a factor of 
200, and they populate the steep, faint end of the 
infrared luminosity function, with Lm < 10 Lq. 
However, IR luminous galaxies, including SMGs 
and ones from our sample, are probably more lu- 
minous by a factor of > 100. Therefore, the inte- 
grated luminosity density from SMGs and Spitzer 
selected ULIRGs at z ~ 2 is roughly comparable 
to that from UV selected galaxies. 

6. Summary 

We have used the IRS on Spitzer to observe 
a sample of 52, 24 /im brighter than lmJy galax- 
ies in the XFLS. Of the 52 targets in our sam- 
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Fig. 10. — Our sample compared with a number of 
other well-known samples including: UV-selected 
(UGRs), near-IR selected (BzKs, DRGs), far-IR- 
selected (SMGs), and near-IR-selected (DRGs). 
The other survey data compilation is based on 
Fig. 11 in Reddy et al. (2006). To derive the 
values for our sample we assumed a flat UV-slope, 
and log(ZiR,) = log £5.8 + 1.4. The solid line shows 
the best-fit relation for UGR galaxies (Reddy et 
al. 2006) , while the dashed line shows the best-fit 
line for the z = sample. 

pie, 47 (90%) have measurable redshifts based 
on the IRS spectra, with 35 at 1.5 < z < 3.2, 
and 12 at 0.62 < z < 1.3. The observed spec- 
tra fall crudely into three categories, one third 
with strong PAH emission, one third with only 
deep silicate absorption and the remaining third 
of the sample with strong mid-IR continuum plus 
weak emission and/or absorption features. The 
sources in our sample are very luminous, with 
estimated 5.8 /jm (i^L„(5.8/im)) continuum lumi- 
nosity > 1O 12 L0, implying Lir on the order of 
> 10 13 L©. Our sources, by selection, have very 
high infrared to the rest-frame UV luminosity ra- 
tios (L m /L ° - 100 - 1000), and are, thus 
v x/ 1600A n ' 

heavily dust obscured in the UV. At the mid-IR 
10 /im, ^60% of our sample have optical depth 
79.8/OTi greater than 1, implying high dust obscu- 
ration along line-of-sight. Over all, our sources 
are among the most luminous, and on average, 
the most dust obscured at z ~ 2. Our survey re- 
vealed an infrared luminous population at z ~ 2, 



which has very little overlap with SMGs and other 
optical/near-IR selected galaxies. 

For the third of the sample with PAH emission, 
we estimate the total energy in the PAH features 
(continuum subtracted) is roughly 20-70% of the 
5.8 /mi continuum luminosity (^L„(5.8/um)), cor- 
responding to ~ 1 - 3% of the total infrared lumi- 
nosity. This implies that PAH features are impor- 
tant components of mid-IR SEDs at high redshift. 
On average, our PAH emitting ULIRGs are about 
a factor of (5-10) more luminous than starburst 
ULIRGs at z ~ 0. The higher luminosity from 
these systems could be due to the combination of 
very efficient star formation and higher fraction of 
gas-to-mass ratio at z ~ 2. Our study has revealed 
a population of extremely IR luminous galaxies at 
z r-j 2, whose physical properties are diverse, with 
a third of the sample powered by dusty star forma- 
tion, one third of sample being deeply embedded 
systems, whose dust heating sources are heavily 
obscured by dust and their energetic nature can 
not be determined unambiguously by the present 
dataset alone. Finally, the remainder third of the 
sample are probably AGN powered systems. The 
enormous energy output and high fraction of dust 
content make these sources prime candidates for 
long-sought type II quasars. 
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MIPS180 


SST24 


J171543 


54+583531 


2 


1683.941 


31.977 





896 




8x2x122 


8x2x122 


MIPS8135 


SST24 


J171455 


71+600822 


6 


1769.055 


103.991 





294 




8x2x122 


8x2x122 


MIPS 15840 


SST24 


J171922 


40+600500 


4 


1820.508 


188.874 





462 


1x2x242 


6x2x122 


7x2x122 


MIPS133 


SST24 


J171433 


68+592119 


3 


1989.473 


150.865 


4 


173 


1x2x242 


2x2x122 


2x2x122 


MIPS22303 


SST24 


J171848 


80+585115 


1 


2029.962 


95.502 





395 




2x2x122 


3x2x122 


MIPS22277 


SST24 


J171826 


67+584242 


1 


2298.740 


227.975 


1 


856 


1x2x242 


2x2x122 


3x2x122 


MIPS110 


SST24 


J171215 


44+585227 


9 


2300.304 


261.559 


1 


890 


1x2x242 


2x2x122 


3x2x122 


MIPS8034 


SST24 


J171210 


28+601858 


1 


2909.007 


193.751 


4 


251 


1x2x242 


2x2x122 


3x2x122 


MIPS78 


SST24 


J171538 


18+592540 


1 


2972.799 


252.988 





186 


1x2x242 


2x2x122 


3x2x122 


MIPS22204 


SST24 


J171844 


38+592000 


5 


4100.928 


428.290 


2 


996 


1x2x242 


2x2x122 


3x2x122 


MIPS42 


SST24 


J171758 


44+592816 


8 


4711.865 


632.269 





186 


1x2x242 


2x2x122 


3x2x122 



a For 8 /J.m none- detected sources, we list la flux limits, taken from the XFLS IRAC catalog bv lLacy et"al| ||2005fl . 

b SL and LL stand for Short-Low and Long-Low module of the IRS. The total exposure time is equal to No. of 
repeats X 2 (No. of dither position) X individual exposure time. 
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Table 2 

The Measured Redshifts and Types for Our Sample 



MID 


%vis i If 


log £5.8 




Type 


features 


A/TTPQ99/10/1 
IVllJr DZZ4U4 


nfii ±n oa 

U.Dl zt U.UO 


1 A Q9-UO 07 
lU.oZztU.U ( 


a 


1 
1 


noil q 10 Q 

y.o,ii.o,iz.o 


T\/TTPQ99^/1 
IvlLr oZZ004 


n S9 -U n nc^ 
U.oZ zt u.uo 


1 n c;o-uo 00 
lU.0Uztu.Uy 


a 


1 
1 


0.0, y. 0,11. o,iz.o 


A,TTT3CQOn7 

MlroozU < 


o s/i j- n aq 
U.o4 zt U.Uo 


1 n 71 J-O 07 
1U. ( IztU.U ( 


a 


1 


y.o,ii.o,iz.o 


iv±_l_t qzzuuu 


n of. _i_ n no 
U.OU ZIZ u.uo 


If) C.Q-UO 10 
iu. oyzcu. ±\j 


cl 


1 

_L 




IVllJr oZoo 


U.y4 ZIZ U.U4 


1 n S^4-0 OS 
lU.oOztU. Uo 


a 


1 
1 


QS 11 Q IO O 


a/ttpqi aoqo 

IVlUrolOUoU 


n os -u n H/i 
u.yo zt u.u4 


1 n QQ-l-0 07 

lu.yoztu.u * 


a 


1 
1 


0811 ^(198 


TVTTPQS1 OA 

lVllJroolo4 


n Qo -u n 09 
u.yy zt u.uz 


1 n Q/i-i-0 oa 
1u.y4ztu.uo 


a 


1 
1 


noil 190 u t 
y. 0, 11. o,iz. 0,10.0 


A/TTPm ^Q9S 
iviiJr oioyzo 


i c;o _i_ n 09 
1 .OZ Zt u.uz 


1 1 /L/L4-0 0^ 
1 1 .44ZtU.UO 


a 


1 
1 


fl 9 7 78 fi 11 ^198 


IVllJr uZZOOl 


i 70 _|_ n 09 
i . ( o zt u.uz 


1 1 .ouztu. 10 


a 


1 
1 


fi 9 7 7 S fi 11 q 19 S 


l\/TTPt;99977 
iviiJr dzzz * * 


i 77 _i_ n 0*3 
1 . f f Zt u.uo 


1 9 0°.4-0 OS 
IZ.UOztU.Uo 


a 


1 
1 


P. 9 7 7 Q 8 
d.z, / . ( ,y .0 


A/TTP^S/LCH 
IVllJr oo4yo 


i on _|- n 01 
l.oU zt U.U1 


I 1 QQ-UO 9°. 

I I .OoZtU.ZO 


a 


1 
1 


77Q811 Q 19 7 

/ . ( ,y .0,1 i.o,iz. ( 


A/TTPQ99/1 S9 
IvllJro ZZ4oZ 


i oa _i_ n 09 

1 .o4 zt U.UZ 


1171 -1-0 OS 
11. / IztU.Uo 


a 


I 
1 


7 7 8 fi Q 8 1 1 1 

/ . * ,o.u,y. 0,11.0 


MTPQ9SQ 
IVllJr OZoj 


i of\ _i_ n 01 

l.oO zt U.U1 


11/1 7-1-0 9Q 
11.4 / ztU. Zo 


a 


1 
1 


fi9778KQ811 

D.z,/./,o.o,y. 0,11.0 


A/TTPQ99^QO 
IVlUroZZOoU 


i _i_ n 01 
i.yo zt u.ui 


1 1 A S-l-0 9^. 
1 1 .4oztU.Z0 


a 


I 
1 


R9778fiQ811 1 
D.Z,/./,o.0,y. 0,11.0 


A/TTPCM 90 

ivnio4zy 


o no -i- n 09 
z.uy zt u.uz 


1 1 fifi_i_n oq 
1 i.ooztu.uy 


a 


1 
1 


ArTT 7 7 8 Q 8 
Aril, / . < ,0.0, y .0 


l\ ; fTpci £1 A A 


o i q _i_ n 09 
Z. 10 Zt u.uz 


1 1 c;Q-|_o 9S 
1 1 .oyztu. zo 


a 


1 

1 


fi9778fiQ811 "? 

o.z, / . ( ,0.0, y. 0,11.0 


MTPQ^OA 
IVllJr oOUO 


9 co_l n no 
Z.OZ zt u.uz 


I 1 Q9_|_n 90 

I I .oZztU.ZU 


a 


1 
1 


fi 9 7 7 Q 8 
O.z, / . ( ,y.o 


IVllJrolOooU 


1 CZA -U O O^ 

1.04 zt U.UO 


1 1 O/IJ-O OQ 

H.y4ztU.Uo 


a 


1 K 

1.0 


1 7 O Q 

( . ( ,y.o 


l\/rTPQ/l A /I 
1V111^o4D4 


1 OK -U O 07 
I.00 zt U.U / 


1 1 KQ-UA OA 
1 l.OoztU.UO 


a 


1 K 
1.0 


7 7 


A/TTPQ1 fil 
IVlLrolOllo 


1 OQ -U n 07 
l.yo zt U.U ( 


1 1 77_|_n OS 
1 1. (t ztU.Uo 


a 


1 K 
1.0 


7 7 Q 
1.1,9.0 


MTPS49 

1V11JT UTrZi 


i i n 07 
i . yo ziz u.u* 


1 ci _la no 

J_ Z . O J- ZIZ U.UZ 


a 


1 ^ 

1 . 


7 7 Q 8 
1 . 1 , 9.0 


A/TTPQ1 A1 99 
IvlUrolOlZZ 


i Q7 _i_ n 0^ 

i.y i zt u.uo 


11 01 _i_n r\p. 
1 l.olztU.UO 


a 


1 ^ 
1.0 


7 7 Q 8 1 1 1 

1 . t ,y. 0,11.0 


a/ttpqi aoso 

IVllJrQlOUoU 


9 0/1-1-0 OA 
Z.U4 zt U.UO 


1 1 QQ_i_n n/i 
1 l.yoztU.U4 


a 


1 ^ 
1.0 


7 7 Q 8 

1 . / ,y.o 


A/TTPQ9990/1 
IVllJr DZZZU4 


9 OS -1- OQ 
Z.Uo zt U.Uo 


1 9 £04-0 09 
IZ.OUztU.UZ 


a 


1 ^ 
1.0 


7 7 Q 8 1 1 1 

1 . / ,y. 0,11.0 


A/TTPm AO^Q 
iviiiuiouoy 


9 /L°. -1- 07 
Z.40 Zt U.U ( 


I 1 AQ4-0 0^ 

I I .oyztu. uo 


a 


1 ^ 

1 .0 


fi 9 7 7 Q 8 
o.z, 1 . / ,y.o 


iv/ripci so 
IVlll^oloU 


9 /1 7 _|_ n 0/1 
Z.4 i zt U.U4 


1 9 Qo_|_n no 
IZ.oZztU.Uo 


a 


1 ^ 
1.0 


fi 9 7 7 Q 8 
O.z, / . / ,y.o 


A/TTPSSQ97 


9 AS + OA 

Zi.4:0 ZIZ U.UU 


1 O 07-UO 07 

1 .U ( Z1ZU.U ( 


a 


1 ^ 

1 . 


fi 9 7 7 Q 8 
u . ^ , 1 . ( , y .0 


MIPS8196 


2.6 ±0.1 


12.46±0.02 


a 


1.5 


7.7,9.8 


MIPS8242 


2.45 ±0.04 


12.07zb0.04 


a 


1.5 


6.2,7.7,8.6,9.8 


MIPS22699 


2.59 ±0.04 


12.12zt0.04 


a 


2 


6.2,7.7,8.6,9.8 


MIPS78 


2.65 ±0.1 


12.66zb0.02 


a 


2 


9.8 


MIPS22303 


2.34 ±0.14 


12.32±0.06 


a 


2 


7.7,9.8 


MIPS8245 


2.7 ±0.1 


12.23±0.08 


a 


2 


9.8 


MIPS22558 


3.2 ±0.1 


12.49±0.02 


a 


2 


7.7,9.8 


MIPS8342 


1.57±0.09 


11.50zt0.08 


a 


3 


6.2,7.7,11.3,12.8 


MIPS22661 


1.75 ±0.03 


11.44zb0.07 


a 


3 


7.7,11.3 


MIPS16095 


1.81 ±0.05 


11.69zb0.05 


a 


3 


9.8,11.3,12.8 


MIPS227 


1.87 ±0.05 


11.68zb0.06 


a 


3 


7.7,9.8,11.3,12.8 


MIPS15958 


1.97 ±0.05 


11.95zt0.05 


a 


3 


7.7,9.8 


MIPS15949 


2.15 ±0.03 


11.94zb0.04 


a 


3 


7.7,11.3 
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Table 2 — Continued 



MID 




log £ 5 .8 


Qz 


Type 


features 


MIPS15840 


2.3 ±0.1 


12.28±0.03 


a 


3 


weak 9.8 


MIPS15977 


1.85 ±0.07 


11.82±0.05 


b 


3 


weak 6.2,7.7,11.3 


MIPS110 


1.0 ±0.1 


11.21±0.07 


b 


3 


10 /im 


MIPS133 


1.0 ±0.2 


10.97±0.05 


b 


3 


10 


MIPS22467 


0.7 ±0.2 


10.57±0.05 


b 


3 


9.8 


MIPS8268 


0.8 ±0.2 


10.82±0.05 


b 


3 


10 /im 


MIPS8034 


0.7 ±0.2 


11.21±0.08 


b 


3 


10 


MIPS279 


999 a 


11857664 


a 


4 




MIPS22582 


999 


11870208 


a 


4 




MIPS15929 


999 


11865088 


a 


4 




MIPS16006 


999 


11856896 


a 


4 




MIPS8135 


999 


11861248 


a 


4 





a The value 999 means that no redshifts were derived based on the IRS 
spectra. 



Table 3 

The Keck Spectroscopy of the IRS Sample 



MID 


z keck 


Z/KS 


Type(IRS) 




Type(Kcck) 


NIRSPEC Data 


MIPS22530 


1.9511 


1.96 


strong PAH (t = 


= 1) 


Ha 


MIPS8327 


2.4426 


2.48 


highly obscured (i 


= 1.5) 


Ha 


MIPS506 


2.4691 


2.52 


strong PAH (t = 


= 1) 


Ha 


MIPS16059 


2.3252 


2.43 


highly obscured (t 


= 1.5) 


Ha 


MIPS22204 


1.9707 


2.08 


highly obscured (t 


= 1.5) 


Ha 


MIPS16080 


2.0077 


2.04 


highly obscured (i 


= 1.5) 


Sy2, [OIII]5007, Ha 


LRIS Data 


MIPS8207 


0.8343 


0.84 


strong PAH (t = 


= 1) 


Sy2, [OII]3727,[NeIII]3868 


MIPS133 


0.9060 


1.0 


weak features (t 


= 3) 


Sy2, [OII]3727,[NeIII]3868,[OIII]5007 


MIPS283 


0.9374 


0.94 


strong PAH (t = 


= 1) 


Sy2, [OII]3727,[NeIII]3868 


MIPS227 


1.6370 


1.87 


weak features (f. 


= 3) 


AGN, Lya,CIV], Hell, CHI], MgH 


MIPS16030 


0.9865 


0.98 


strong PAH (t = 


= 1) 


star forming, [Oil] 3727 


MIPS16080 


2.011 


1.95 


highly obscured (t 


= 1.5) 


QSO, broad CII]2326 


DEIMOS Data 


MIPS22404 


0.6101 


0.61 


strong PAH (t = 


= 1) 


star forming, [Oil] 3727 


MIPS22600 


0.8504 


0.86 


strong PAH (t = 


= 1) 


star forming, [Oil] 3727 


MIPS22554 


0.8242 


0.82 


strong PAH (t = 


= 1) 


star forming, [Oil] 3727 


MIPS279 


1.2319 


1.3 


weak features (f. 


= 3) 


AGN 


MIPS110 


1.0505 


1.0 


weak features (t 


= 3) 


AGN 
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